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ABSTRACT 

Observations of the damped Lya systems provide direct measurements on the chemical enrichment 
history of neutral gas in the early universe. In this Letter, we present new measurements for four damped 
Lya systems at high redshift. Combining these data with [Fe/H] values culled from the literature, we 
investigate the metallicity evolution of the universe from z w 1.5 — 4.5. Contrary to our expectations 
and the predictions of essentially every chemical evolution model, the A(HI)-weighted mean [Fe/H] 
metallicity exhibits minimal evolution over this epoch. For the individual systems, we report tentative 
evidence for an evolution in the unweighted [Fe/H] mean and the scatter in [Fe/H] with the higher 
redshift systems showing lower scatter and lower typical [Fe/H] values. We also note that no damped 
Lya system has [Fe/H] < — 2.7 dex. Finally, we discuss the potential impact of small number statistics 
and dust on our conclusions and consider the implications of these results on chemical evolution in the 
early universe. 

Subject headings: galaxies: abundances — galaxies: chemical evolution — quasars : absorption lines 



1. INTRODUCTION 

The damped Lya systems, neutral hydrogen gas lay- 
ers identified in the absorption line spectra of background 
quasars, dominate the neutral hydrogen content of the 
Universe at all epochs. At high redshift, these systems are 
widely accepted as the progenitors of present-day galaxies 



metallicity of the damped Lya systems from z ps 2 — 4.5 
which places tight constraints on chemical evolution mod- 
els (e.g. Pei et al. 1999), as well as a valu able consistency 
check on SFR observations (Pettini 1999). 

Over the past decade, several groups have surveyed the 
metallicity of the damped Lya systems from z ~ 1—4 
dPettini et al. 1994 |Pettini et al. 1997] [Pettini 1999| , |Lu ct] 



10 — cm , imply over 



densities Sp/p ^> 100, i.e., these are virialized systems at 
high redshift; (ii) they contain the majority of neutral gas 
in the early Universe and are therefore the reservoirs for 
galaxy formation; (iii) their gas density fi gos at redshift 
z w 2 — 3 is consistent with the mass density of stars to- 
day (Wolfe et al. 1995). While the physical nature of the 
damped Lya systems is still controversial ( Prochaska fc 



for th e following rcaaonu: (i) their very large HI column ' I d . r — \ ' u - mnnh rJ S— r — rr — r 1 • — i 

j ... at/ttt\ ^ at o m20 -2 ■ i al. 1996 , Prochaska & Wolfe 1999 . To date, the chemical 
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Wolfe 1997j lHachnelt ct al. 1998| , |Maller et al. 1999| , [Le 



Brim ct al. 1996), by studying the chemical abundances 



of the damped Lya system one directly traces the chem- 
ical enrichment history of the Universe at high redshift. 
Observing damped Lya systems is equivalent to poking 
sightlines through the ISM of protogalaxies. Because these 
observations are biased to HI cross-section and the HI gas 
mass of a system is proportional to J a ■ N, one can mea- 
sure global properties of the universe simply by weighting 
the measurement from each damped system by A(HI). At 
the same time, the observations afford an efficient means 
for examining the characteristics of individual protogalax- 
ies in the early Universe. In this Letter, we examine the 

1 Visiting Astronomer, W.M. Keck Telescope. The Keck Obser- 
vatory is a joint facility of the University of California and the Cali- 
fornia Institute of Technology. 



abundances of over 40 systems have been measured, the 
majority with z — 1.5 — 3 where the identification and 
follow-up observations of damped Lya systems is most 
efficient. These studies argue that at z w 2, the mean 
metallicity of the damped systems is approximately 1/10 
- 1/30 solar metallicity ([Zn/H] « -1.1, [Fe/H] « -1.5) 
with a large scatter from nearly solar to less than 1/100 
solar metallicity. At very high redshift (z > 3), the picture 
is far less certain. Focusing on a sample of seven z > 3 
damped Lya systems, Lu et al. (1996,1997) noted that the 
metallicity of these systems is significantly lower than the 
z < 3 observations. In turn, the authors argued that z ss 3 
marked the epoch where significant star formation begins, 
a claim with important implications for the processes of 
galaxy formation. 

In this Letter we present new measurements on the 
metallicity of four damped Lya systems (including three 
at z > 3.5) and together with the data from Prochaska 
& Wolfe (1999) double the sample of z > 3 systems. 
The new full sample - including the systems from Lu et 
al. (1996,1997) - reveals evidence for little change in the 
A(HI)-weighted mean metallicity of the neutral Universe 
from z w 1.5—4.5, contrary to the predictions of essentially 
every chemical evolution model. On the other hand, we 



2 



METALLICITY EVOLUTION IN THE EARLY UNIVERSE 



find tentative evidence for an evolution in the unweighted 
mean and scatter of [Fe/H] for individual damped Lya 
systems. Finally, we comment on the robustness of these 
results (particularly in the light of small number statis- 
tics and dust), speculate on the implications for chemical 
enrichment, and discuss the prospects for future advances. 

2. OBSERVATIONS AND ANALYSIS 

To determine the metallicity of a damped Lya system, 
one must accurately measure the neutral hydrogen col- 
umn density iV(HI) and a metallicity indicator, typically 
either Zn or Fe. In stellar population studies of the Galaxy 
one traditionally uses Fe as the metallicity indicator, pri- 
maril y as a matter of convenience. As we are studying 



Table 1 

new metallicity measurements 



QSO 


Z-abs 


AT (HI) 


Af(Fc+) 


[Fe/H] 


BRI0952-0115 


4.024 


20.50 ±0.1 


14.054 ± 0.07 a 


-1.95 


BRI1108-0747 


3.608 


20.50 ±0.1 


13.860 ± 0.03 b 


-2.14 


Q1223+1753 


2.466 


21.50 ±0.1 


15.279 ± 0.03 c 


-1.72 


PSS1443+2724 


4.224 


20.80 ±0.1 


15.325 ±0.10 c 


-0.98 



a Average of Fell 1144 (logs/ = 0.105) and Fell 1608 
6 FeII 1608 
c FeII 1611 



were derived primarily from the Fell 1608 and/or Fell 1611 



gas-phase abundances, however, we must account for the 



transitions with the ap parent optical depth method (Sav 



possible depletion of Fe onto dust grains or instead choose 
an element like Zn which is minimally affected by deple- 
tion. Unfortunately, there are both theoretical and ob- 
serva tional disadvantages to using Zn as the metallicity 



indicator. Theoretically, Zn has a very uncertain chemical 
origin. It is referred to as a n Fe peak element b ecause it 
traces Fe in Galactic stars (Sneden et al. 1991), yet the 
leading theory on the production of Zn propo ses it forms 
in the n eutrino-driven winds of Type II SN (Hoffman et 
al. 1996| ) . Furtherm ore, recent me asurements of |Zn/Fc| 
i n metal-poor stars (|Johnson 199~ ) and thick disk stars 



( Prochaska et al. 206o| ) suggest Zn/Fe is enhanced rela- 
tive to the Sun by +0.1 to +0.3 dex, perhaps consistent 
with a Type II origin. Observationally there are compli- 
cations with measuring Zn in the damped Lya systems, 
where one must rely on two weak Znll transitions with 
Xrest ~ 2000A. The transitions are so weak that even at 
high resolution and high S/N, Zn can only be detected in 
damped systems when log[V(HI)] + [Zn/H] > 19.0 (e.g. 
[Zn/H] > -1.3 for systems with AT(HI) w N thresh ). Most 
important to this study, however, the large rest wavelength 
of the Znll transitions prevents one from readily measuring 
Zn in z > 3 damped Lya systems as it is difficult to make 
sensitive observations at A w 8000A with current high res- 
olution spectrographs. In fact, at the time of publication 
we are not aware of a single accurate Zn measurement for 
any z > 3 damped Lya system. Therefore, we will focus 
on Fe in this Letter, which has two singly ionized tran- 
sitions at X re st ~ 1600A with a complement of /-values 
ideal for measuring the abundance of Fe in damped sys- 
tems. We restrict the analysis to Fe measuremen ts made 
with HIRES on the Keck I telescope ( |Vogt 1992j ), specif- 
ically the systems observed by Prochaska & Wolfe (1999) 
and Lu et al. (1996,1997) and the additional systems in- 
troduced here. In addition to providing a homogeneous 
data set which has been reduced and analyzed with the 
same techniques, these observations account for nearly ev- 
ery damped Lya system with an accurate Fe abundance 
at z > 1.5 and every system with z > 3. 

Table [l] summarizes the new V(Fe + ) measurements de- 
rived from observations acquired by the authors in Febru- 
ary 1998 and March 1999 with HIRES on the Keck I 10m 
telescope. The data was reduced with the makee software 
package developed by T. Barlow and the column densities 



age and Scmbach 1991). We adopt the oscillator strengths 
from Cardclli & Savage (1995) noting that our conclu- 
sions on the evolution of [Fe/H] in the damped Lya sys- 
tems are not sensitive to their values. The iV(HI) val- 



ues for these systems arc taken from the liter ature (Wolfe 
et al. 1995| , [Btorric-Lombardi fc Wolfe 2000D and are the 
dominant source of error in the [Fe/HJ values. Finally, 
we evaluate [Fe/H] assuming the meteoritic Fe abundance 
(e(Fe)=7.50; Grevesse & Sauval 1999) without adopting 
any ionization corrections which is an excellent assumption 
for all but possibly the low est iV(HI) damped Lya systems 
(Prochaska & Wolfe 1996). Together with the published 
measurements of Prochaska & Wolfe (1999) and Lu et al. 
(1996,1997) the total [Fe/H] sample is 37 systems, 15 with 
z > 3. The systems were chosen independent of any prior 
metallicity measurements; the only possible biases are due 
to the magnitude limited selection of the quasars (e.g. Fall 
& Pei 1993) which will be discussed in the following sec- 
tion. 



3. RESULTS AND DISCUSSION 

Figure [l] plots the 37 [Fe/H] values versus z a t, s for the 
Wolfe & Prochaska sample (dark squares) and (light stars) 
the sample of damped Lya systems observed by Lu et al. 
(1996,1997). To explore evolution in the metallicity of the 
damped Lya systems, we consider three moments of the 
metallicity data in two redshift intervals, zi ow = [1.5,3] 
and Zhigh = (3,4.5]. These are: (1) the V(HI)-weighted 
mean metallicity of neutral gas, < Z >; (2) the unweighted 
mean metallicity, < [Fe/H] >, of the set of damped Lya 
systems at zi ow and Zhigh] an< i (3) the standard deviation 
of [Fe/H] in these protogalaxies, <r([Fe/H]). 

The first moment represents the global metallicity of all 
neutral gas at a given epoch, £l m etais /^hi- It is evaluated 
by weighting each [Fe/H] measurement by the correspond- 
ing HI column density, < Z > = log[£N(Fe+)/Y,N(m)}- 
log(Fe/H) Q . Computing the mean for the damped Lya 
systems at the two epochs, we find < Z >i ow = —1.532 ± 
0.036 and < Z >u g h= -1.634 ± 0.049. The errors on 
the < Z > values reflect only the statistical uncertainty 
in measuring iV(Fe + ) and AT(HI) and were derived with 
standard error propagation techniques. Below we estimate 
the uncertainty due to small number statistics. Compar- 
ing the < Z > values, we note that they favor no signifi- 
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Fig. 1. — Thirty-nine [Fe/H], z a bs pairs for the damped Lya systems observed by Prochaska & Wolfe (1999, this paper; 
squares) and Lu et al. (1996,1997; stars) with HIRES on the Keck I telescope. The open circles correspond to the iV(HI)- 
weighted mean metallicity for the systems at z a b s = [1.5,3] and z a i> s — (3,4.5]. Note that the difference in these means is 
small, k, 0.1 dex. Also observe that the scatter in the [Fe/H] values appears to increase at lower redshift. 



cant evolution in the mean metallicity of neutral gas from 
z = 1.5 — 4.5. If we include the tentative result from Pet- 
tini et al. (2000) that the Zn mean metallicity does not 
change from z w 1 — 3, then one concludes there is no evi- 
dence for significant metallicity evolution from z — 1 — 4.5, 
an interval spanning more than 3 Gyr. The other two 

n 

moments, the unweighted mean < [Fe/H] >= i^[Fe/H] 
and the scatter <r([Fe/H]), are more sensitive to the chem- 
ical enrichment history within individual protogalaxies as 
each damped system is given equal weight. For the two 
intervals we find that the mean logarithmic abundance, 
< [Fe/H] > low = -1.61 and < [Fe/H] > high = -1.83. 
Meanwhile, the scatter in [Fe/H] is cr([Fe/H]) = 0.50 and 
er([Fe/H]) = 0.35 for the zi ow and Zhigh samples respec- 
tively. Performing the Student's t-test and the F-test 
on the two moments, we find that the < [Fc/H] > and 
<r([Fe/H]) statistics for the two epochs are inconsistent at 
the 90% and 80% c.l. Therefore, there is tentative evi- 
dence for chemical evolution in the individual damped Lya 
systems with the z < 3 sample showing a higher typical 
metallicity and a larger scatter in [Fe/H] from system to 
system. 

To address the robustness of these results, one must con- 
sider several issues. First, because < Z > is dominated by 



the s ystems with the largest iV(III) values, this mean is 



robus t only in do far ao the total iV(HI), Rl T - ]P AT (HI), 

n 

well exceeds that of a single damped Lya system. Fig- 



ure g plots the [Fe/H], AT(HI) pairs for all 37 systems 
where the dark circles are members of the zi ow sample 
and the light triangles those of the Zhigh sample. Note 
that there are four systems with A(HI) > 10 21 cm" 2 at 
z < 3 but only a single system in the Zhigh sample. For 
the ziow sample, Hly = 10 22 30 cm -2 , which is a fac- 
tor of 4 larger than the largest iV(HI) measured for any 
damped Lya system (Q0458-0203; AT(HI) = 10 217 ) and 
10 times greater than most of the known damped Lya 
systems. As such, we consider the mean derived from 
the zi ow sample to be reasonably robust. The primary 
potential pitfall is if the optical surveys have systemati- 
cally missed damped systems with iV(HI) > 10 22 cm -2 , a 
possibility if dust obscuration is significant (discussed fur- 
ther below). The situation is far more uncertain for the 



z > 3 sample where Hly 



10 



21.84 



cm , comparable to 



the A(HI) of the Q0458— 0203 system from the Zi ow sample 
and only three times greater than the largest iV(HI) sys- 
tem in the Zhigh sample. While the most recent surveys 
suggest there are very few z > 3 damped systems with 
AT (HI) > 10 21 cuT 2 Storrie-Lombardi & Wolfe (2000), we 
caution that the mean we have derived for the Zhigh sam- 
ple is a tentative result. For example, the system toward 
Q0000— 2619 has significant bearing on < Z >h%gh and its 



Fe abundance has been diffi cult to determine (Prochaska 



& Wolfe 199E, Lu et al. 1996). Ironically, removing it from 
the z^gh sample would actually increase < Z >high into 
exact agreement with < Z >i ow because we have adopted 
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Fig. 2. — Thirty-nine [Fe/H], V(HI) pairs for the damped Lya systems in the full sample. The circles correspond to 
z a bs < 3 systems and the triangles are for z a f, s > 3. While the systems with AT(HI) < 10 21 cm~ 2 show a large scatter in 
[Fe/H], the large AT (HI) systems all have [Fe/H] w —1.8. 



[Fe/H] = -1.77 based on the Fell 1611 profile from this 
system. Meanwhile, lowering [Fe/H] by 0.6 dex to estab- 
lish consistency with the Ni and Cr abundances, would 
decrease < Z >high by 0.1 dex. In short, while we have 
confidence in the < Z >i ow value, we caution the reader 
that small number statistics are still important in evaluat- 
ing the < Z >high- One can estimate the uncertainty asso- 
ciated with the small number statistics of the two samples 
by performing a bootstrap statistical analysis. For each 
sample, we independently calculated < Z > 500 times by 
randomly drawing n objects (n is the number of damped 
systems in a given redshift interval) from each interval. In 
turn, we can estimate the effects of cosmic variance on our 
results by calculating the standard deviation of the two 
bootstrap < Z > distributions: o-f- z> = 0.088 dex and 



<z> 

high 



0.155 dex. As one would expect, the < Z > 



value, which is based on only 15 systems, is considerably 
less certain than the < Z >i ow measurement. The differ- 
ence in the a <z> values stresses the outstanding need for 
future observational programs to focus on z > 3 damped 
systems. 



Any study on the chemical abundances of the damped 
Lya systems must assess the potential effects of dust. 
With respect to this analysis, where we have taken Fe 
as the metallicity indicator, there are two important as- 
pects to consider: (1) if we need to correct the observed 
[Fe/H] values by some factor to obtain the true metal- 
licity of each system, does the mean correction evolve in 
time and/or differ from system to system at the same 



epoch?; and (2) dust obscuration could remove damped 
Lya systems from the magnitude limited samples which 
would significantly alter the conclusions (e.g. metal-rich, 
high AT (HI) systems). With respect to the first concern, 
we can estimate the maximum dust correction to [Fe/H] 
via the measured Zn/Fe ratio. Again, Zn is essentially 
undepleted in the gas-phase so that [Zn/H] = [Fc/H] + 
[Zn/Fe] may be more representative of the true metallicity 
in the damped Lya systems. This practice is limited, how- 
ever, by the fact that Zn may be produced in Type II SN 
( Hoffman et al. 1996 ) such that super-solar Zn/Fe ratios 
would be representative of nucleosynthesis, not dust deple- 
tion. Furthermore, recent results on the [Zn/Fe] ratio mea- 
sured in Galactic stars shows that [Zn/Fe] > +0.2 dex in 
very metal-poor stars ([Fe/H] < —2.5; Johnson 1999) and 
even exhibits super-solar values (average [Zn/Fe] « +0.13 
in 10 stars; Prochaska et al. 2000) in thick disk stars 
with [Fe/H] > -1. Therefore, while the typical [Zn/Fe] 
value in the dampe d Lya systems is +0-4 dex with rela- 
tively small scatter (Pettini et al. 1997, Prochaska fc Wolfe] 



1999), it is unclear what fraction is due to dust depletion. 



Nonetheless, if we take [Zn/H] as the true metallicity in- 
dicator, < Z > and the unweighted mean are enhanced by 
w 0.4 dex, but there is very little change in the observed 
scatter. The potential effects of dust depletion on the sta- 
tistical moments for the z^igh sample are more speculative 
as there is no accurate Zn determination for any z > 3 
damped Lya system. To estimate the depletion level, we 
can compare the relative abundance patterns (in particular 
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the Si/Fe ratio) of these systems with the zi ow sample. In 
the few cases where Si/Fe has been measured in the z > 3 
systems one finds [Si/Fe] « +0.3 dex, nearly identical to 
the typical value of the z < 3 sample. While the similarity 
of a metal ratio like Si/Fe does not require similar dust 
depletion levels, the z > 3 [Si/Fe] values do imply deple- 
tion levels of at least 0.3 dex. Therefore, unless one takes 
the unlikely stance that the z > 3 systems are significantly 
more depleted than the zi ow sample, we expect very min- 
imal evolution in the depletion levels of the damped sys- 
tems and no significant impact on any of our conclusions. 
The effects of biasing due to dust obscuration are more 
difficult to address. Note in Figure || the absence of any 

AT(HI) > 10 21 cm- 2 systems with [Fe/H] 1. While 

this may be due to small number statistics or that very 
few regions exist in the early universe with large AT(HI) 
and [Fe/H] ^ —1, the trend could also be explained by 
dust obscuration. Fall & Pei (1993) have presented an ex- 
cellent framework for addressing the effects of dust deple- 
tion on damped Lya statistics. Their calculations indicate 
that if the logarithmic scatter in the dust-to-gas ratio k is 
small (less than 1 dex) , then only a small correction to the 
mean optical depth and in turn to < Z > is required (Fall 
1999). For a constant dust-to-metals ratio - implied by the 
nearly constant [Zn/Fe] values - the logarithmic scatter in 
k pa er([Fe/H]) and we have shown <r([Fe/H]) < 0.5 for the 
two samples. Therefore, we expect dust obscuration to 
have a minimal effect (< 0.2 dex) on our results. 

The results in this Letter on the evolution of the metal- 
licity of neutral gas in the Universe and individual pro- 
togalaxies present an unexpected picture. A number of 
groups have estimate d the chemical evolution of ne u tral 
gas a t high redshift (Malaney and Chaboycr 1996|, Ed 
~ ' 'alb 



mund f; fc Phillips 1997 , Pei ct al. 199E ) and essentially ev- 
ery treatment predicts a substantial (> 0.5 dex) increase 
in the mean metallicity from z = 4 to z — 2. While a 
0.5 dex evolution is consistent with our results at the 3cr 
level, the current observations favor a very mild evolution 
in < Z >. If future observations lend further support 
for this conclusion, the theoretical models will require sig- 
nificant revision. Of course, these theoretical treatments 
depend sensitively on a number of factors which are un- 
certain: (i) the star formation rate, (ii) the IMF, (iii) the 
mass distribution of protogalaxies, (iv) the loss of metals 
to the IGM, (v) the yield of various elements, etc. There- 
fore, there is considerable theoretical freedom to bring the 
models into agreement with the observed lack of evolution. 
Nonetheless, the Lyman break galaxies offer incontrovert- 
ible evidence t hat significant star formation is taking place 
from z = 3—4 ( Stcidcl ct al. 199§| ) such that the total metal 
content of the Universe must be increasing. Unless these 
metals are enriching only ionized regions (an unlikely sce- 
nario), then to explain the minimal evolution in < Z > 
the total HI content of the Universe must be increasing at 
nearly the same rate as the metal content. It is intriguing 
to note that this is qualitatively consistent with the evolu- 
tion of Qgas observed by Storrie-Lombardi & Wolfe (2000) 
for the damped Lya systems. 

The other statistical moments are sensitive to the chem- 



ical enrichment history within individual galaxies. Com- 
paring the unweighted mean with the weighted mean we 
find that < [Fe/H] > is less than < Z > at both epochs. 
While the difference is not large (pa 0.1 — 0.2 dex), it does 
highlight the fact that many of the AT (HI) < 10 21 cm" 2 
systems exhibit low metallicity In particular, in the 
z > 3 sample only 2 of 15 damped systems show [Fe/H] 
> — 1.5 dex. One possible explanation for the difference is 
systems which have just formed have preferentially low 
AT (HI) and [Fe/H]. The trend is also suggestive of the 
correlation Cen & Ostriker (1999) find between overden- 
sity and metallicity in their numerical simulations. The 
problem remains, however, in explaining why the high- 
est metallicity systems of the zi ow sample also have low 
AT (HI). Lastly, recall that there is tentative support for 
an evolution in both the scatter and < [Fe/H] > with the 
ziow sample yielding larger values. If more recently formed 
systems tend to have lower metallicity, then the evolution 
may easily be explained by a larger mean and scatter in 
the age of the damped Lya systems at z ~ 2. Further- 
more, the systems at z ~ 2 may have larger masses and a 
greater variety of morphologies. 

Finally, we stress that only two systems from the full 
sample have [Fe/H] < —2.5 and the large majority show 
[Fe/H] > -2. As first noted by Lu et al. (1997), there ap- 
pears to be a threshold to the minimum metallicity of the 
damped Lya systems at pa 1/100 solar metallicity There- 
fore, even out to z w 4.5 there is no evidence for damped 
Lya systems with primordial abundance. This places a 
further constraint on chemical evolution models. As we 
probe higher and higher redshift without detecting pri- 
mordial gas, one may be forced toward one of the follow- 
ing conclusions: (i) star formation proceeds rapidly (< 10 
years) to bring the metallicity to 1/100 solar after the for- 
mation of a damped system; (ii) either the damped system 
or its progenitors have been undergoing star formation for 
a lengthy time; and/or (iii) a generation of Population III 
stars has pre-enriched all of the gas. 

The future prospects for improving the z > 3 observa- 
tional sample are excellent. While further progress with 
HIRES and similar instruments is limited by the faintness 
of z > 4 quasars, the new ESI instrument on Keck II will 
be ideal for surveying the AT(HI) and metal content of a 
large (N > 20) sample of very high redshift damped Lya 
systems. We intend to pursue such a program over the 
next few years, taking full advantag e of the ever inc reas- 
ing sample of known z > 4 quasars ( Fan et al. 1999a ) . 



We would like to thank A. McWilliam, E. Gawiser, M. 
Pettini, and M. Fall for insightful discussion and com- 
ments. We thank T. Barlow for providing the HIRES re- 
duction package. We acknowledge the very helpful Keck 
support staff for their efforts in performing these observa- 
tions. J.X.P. acknowledges support from a Carnegie post- 
doctoral fellowship. 



REFERENCES 
Cardelli, J.A. and Savage, B.D. 1995, ApJ, 452, 275 



G 



METALLICITY EVOLUTION IN THE EARLY UNIVERSE 



Cen, R. & Ostriker, J.P. 1999, ApJ, 519, 109L 
Edmunds, M.G. & Phillips, S. 1997, MNRAS, 292, 733 
Fall, S.M. & Pei, Y.C. 1993, ApJ, 402, 479 

Fall, S.M. 1999, private communication . 

Fan, X_, — SDSS collaboration, 1999a, AJ, submitted (astro 

ph/i 903237 ) 

Grevesse, IN. L Sauval, A.J. 1999, a, 347, 348 

Haehnelt, M.G., Steinmetz, M. & Rauch, M. 1998, ApJ, 495, 647 
Hoffman, R.D. et al. 1996, ApJ, 460, 478 
Johnson, J. 1999, PhD thesis, UC Santa Cruz 

Le Brun, V., Bergeron, J., Boisse, P., & Deharveng, J.M. 1997, A&A, 
321 733 

Lu, L.', Sargent, W.L.W., Barlow, T.A., Churchill, C.W., & Vogt, S. 

1996, ApJS, 107, 475 . . 

Lu, L., Sargent, W.L.W., & Barlow, T.A. 1997, astro-ph/971129£ 

Malaney, R.A. and Chaboyer, B. 1996, ApJ, 4627~51 

Mailer, A.H., Somerville, R.S., Prochaska, J.X., & Primack, J.R. 

1999, After the Dark Aq"«- The. Q * h Ajinual October Astrophysics 

Conference, in press ( istro-ph/9812369) 



Pei, Y..C , Fall, S.M., & Hauser, M.G. 191)9, ApJ, submitted, (astro 
ph/i 812182 ) 

Pcttini M., Smith, L. J., Hunstead, R. W., and King, D. L. 1994, 
ApJ, 426, 79 

Pcttini, M., Smith, L.J., King, D.L., & Hunstead, R.W. 1997, ApJ, 
486, 665 

Pcttini, M. 1999, in Lecture Notes in Physics, Proceedings of ESO 

Workshop , od i J i /Walsh and M. Rosa (Berlin: Springer), in press, 

astro-ph/9902173 
Pcttini, M., Ellison, S.L., Steidel, C.C., Shapley, A.E., & Bowen, 

D.V. 2000, ApJ, in press 
Prochaska, J. X. & Wolfe, A. M. 1996, ApJ, 470, 403 
Prochaska, J. X. & Wolfe, A. M. 1997, ApJ, 486, 73 
Prochaska, J. X. & Wolfe, A. M. 1999, ApJS, 121, 369 
Prochaska, J. X. et al. 2000, in preparation 
Savage, B. D. and Sembach, K. R. 1991, ApJ, 379, 245 
Sneden, C, Gratton, R.G., & Crocker, D.A. 1991, A&A, 246, 354 
Steidel, C.C., Adelberger, K., Dickinson, M., Giavalisco, M., Pettini, 

M., & Kellogg, M. 1998, ApJ, 492, 428 
Storrie-Lombardi, L.J. Sz Wolfe, A.M. 2000, ApJ, submitted 
Vogt, S. S. 1992, in ESO Conf. and Workshop Proc. 40, 

High Resolution Spectroscopy with the VLT , ed. M.-H. Ulrich 

(Garching: ESO), p. 223 
Wolfe, A. M., Lanzetta, K. M., Foltz, C. B., and Chaffee, F. H. 1995, 

ApJ, 454, 698 
Wolfe, A.M. & Prochaska, J.X. 1998, ApJ, 494, 15L 



